RANDOM QUANTUM CHANNELS I: GRAPHICAL CALCULUS AND 
THE BELL STATE PHENOMENON 



BENOIT COLLINS AND ION NECHITA 

Abstract. This paper is the first of a series where we study quantum channels from the 
random matrix point of view. We develop a graphical tool that allows us to compute the 
expected moments of the output of a random quantum channel. 

As an application, we study variations of random matrix models introduced by Hayden 
[7], and show that their eigenvalues converge almost surely. 

In particular we obtain for some models sharp improvements on the value of the largest 
eigenvalue, and this is shown in a further work to have new applications to minimal output 
entropy inequalities. 



1. Introduction, motivation & plan 

The theory of random matrices is a field of its own, but whenever it comes to apphcations, 
the driving idea is almost always that although it is very difficult to exhibit matrices having 
specified properties, a suitably chosen random matrix will have very similar properties as the 
original matrix with a high probability. This idea has been used successfully for example in 
operator algebra with the theory of free probability. 

In 2007 for the first time, a similar leitmotiv was used with success by Hayden in [71 and 
Hayden- Winter in [9j to disprove the Renyi entropy additivity conjecture for a wide class of 
parameters p. A proof for the most important case p = 1 was even announced by Hastings 
in [6] with probabilistic arguments of different nature. This is arguably the most important 
conjecture in quantum information theory, and the random matrix models introduced by 
Hayden and their modifications due to Hastings seemed very new from our random matrix 
points of view. 

This paper is therefore an attempt to understand these matrix models with random matrix 
techniques. For this purpose, we introduce a formalism that is very close to that of planar 
algebras of Jones [lO] , and we suggest that the language of quantum gates and planar algebras 
should be considered as very closely related to each other. 

Our paper is organized as follows. In Section [2.11 we recall known facts about integration 
over unitary groups and their large dimension asymptotics. This is nowadays known as Wein- 
garten calculus. In Section [3l we introduce a graphical model to represent (random) matrices 
arising in random quantum calculus. Section |4] gives a theoretical method for computing ex- 
pectations with our graphical model and in the last two sections we give explicit applications 
of these techniques to random quantum channels. More precisely, in Section [5] we investigate 
tensor products of two independent quantum channels and in Section [6] we look at a product 
of a random channel with the channel defined by the conjugate unitary U. 
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Limit theorems presented in this paper are just a sample of what we can be accomphshed 
with the calculus developed in Section HI New results will be given in the forthcoming papers 

[an!- 

2. Background on Weingarten calculus and quantum channels 

2.1. Weingarten calculus. This section contains some basic material on unitary integration 
and Weingarten calculus. A more complete exposition of these matters can be found in [21 [5]. 
We start with the definition of the Weingarten function. 

Definition 2.1. The unitary Weingarten function Wg(n, o") is a function of a dimension 
parameter n and of a permutation a in the symmetric group Sp. It is the inverse of the 
function a ^ n"^" under the convolution for the symmetric group (^a denotes the number of 
cycles of the permutation a). 

Notice that the function a ^ n*'^ is invertible as n is large, as it behaves like n^Se as ra — > oo. 
If n < p the function is not invertible any more, but we can keep this definition by taking 
the pseudo inverse and the theorems below will still hold true (we refer to [5] for historical 
references and further details). We shall use the shorthand notation Wg(o") = Wg(n, a) when 
the dimension parameter n is obvious. 

The function Wg is used to compute integrals with respect to the Haar measure on the 
unitary group (we shall denote by U{n) the unitary group acting on an n-dimensional Hilbert 
space). The first theorem is as follows: 

Theorem 2.2. Let n be a positive integer and i = (ii, . . . , ip), i' = (i[, . . . , i'p), j = (ji, . . . , jp), 
j' = (jj,...,j') be p -tuples of positive integers from {1,2, ... ,n} . Then 



U{n) 



If p ^ p' then 

(2) / u,,,,---u,^,jj^^...w-r^du = o. 

Ju{n) ' ' f f 

Since we shall perform integration over large unitary groups, we are interested in the values 
of the Weingarten function in the limit n — > oo. The following result encloses all the data we 
need for our computations about the asymptotics of the Wg function; see [2] for a proof. 

Theorem 2.3. For a permutation a € Sp, let Cycles{a) denote the set of cycles of a. Then 

(3) Wg{n,a) = {-ir-*^ n Wg(n,c)(l + 0(n~2)) 

c£ Cycles{a) 

and 

(4) Wg(n,(l,...,d)) = (-l)'^-ic,_i H (n-j)-' 



-d+l^j^d-l 



where Ci = ^^^^jj is the i-th Catalan number. 
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A shorthand for this theorem is the introduction of a function Mob on the symmetric group, 
invariant under conjugation and multiphcative over the cycles, satisfying for any permutation 
(J G Sp-. 

(5) Wg(n,f7) =n-(P+l'"l)(Mob(a) +0(n~2)). 

where |cr| = p — i^cr is the length of a, i.e. the minimal number of transpositions that multiply 
to a. We refer to [5j for details about the function Mob. We finish this section by a well 
known lemma which we will use several times towards the end of the paper. This result is 
contained in [12j . 

Lemma 2.4. The function d{a,T) = \a^^T\ is an integer valued distance on Sp. Besides, it 
has the following properties: 

• the diameter of Sp is p — 1; 

• d{-,-) is left and right translation invariant; 

• for three permutations cri,cj2,T E Sp, the quantity d(r, cji) + d(r, CJ2) has the same 
parity as d{ai,a2); 

• the set of geodesic points between the identity permutation id and some permutation 
a & Sp is in bijection with the set of non-crossing partitions smaller than tt, where 
the partition vr encodes the cycle structure of a. Moreover, the preceding bijection 
preserves the lattice structure. 

We end this section by the following definition which generalizes the trace function. For 
some matrices Ai, A2, ■ ■ ■ , Ap £ A4„(C) and some permutation a £ Sp, we define 

Ti^iAi,...,Ap)= Yl TT{Ai,Ai,---Ai^). 

cgCycles(cr) 
c={n i2 ■■■ik) 

We also put Tt^{A) = Tv^{A, A,..., A). 

2.2. Quantum channels. In Quantum Information Theory, a quantum channel is the most 
general transformation of a quantum system. Quantum channels generalize the unitary evo- 
lution of isolated quantum systems to open quantum systems. Mathematically, we recall that 
a quantum channel is a linear completely positive trace preserving map from A^„(C) to it- 
self. The trace preservation condition is necessary since quantum channels should map density 
matrices to density matrices. The complete positivity condition can be stated as 

Vd ^ 1, ^ (Slid ■ ■M.ndi'C) A4„d(C) is a positive map. 

The following two characterizations of quantum channels turn out to be very useful. 

Proposition 2.5. A linear map ^ : A4„(C) A4„(C) is a quantum channel if and only if 
one of the following two equivalent conditions holds. 

(1) (Stinespring dilation) There exists a finite dimensional Hilbert space K. = C^, a 
density matrix Y G Aidi'C) and an unitary operation U € U{nd) such that 

(6) $(X) = = Tr/c [UiX <g) Y)U*] , G 7W„(C). 

(2) (Kraus decomposition) There exists an integer k and matrices Li, . . . , G A^„(C) 
such that 

k 

<^{X) = L^XL*, yx G Mn{C). 

i=l 
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and 

k 

= In • 

1=1 

It can be shown that the dimension of the ancilla space in the Stinespring dilation theorem 
can be chosen d = dim/C = n? and that the state Y can always be considered to be a rank 
one projector. A similar result holds for the number of Kraus operators: one can always find 
a decomposition with k = ii? operators. 

In the final two sections of the paper we study a model of random quantum channels 
originating from the Stinespring dilation formula ([6]). We shall be interested in the spectral 
properties of the elements in the image of such random channels. Quantum channels will be 
the main field of application of the graphical calculus we develop in the next two sections, our 
aim being the treatment of some additivity problems in quantum information theory. 

3. Graphical model 

In this section, we lay out the foundation for the graphical calculus we shall develop later. 
We introduce a graphical formalism for representing tensors and tensor contractions that 
is adapted to quantum information theory. We start at an abstract level, with a purely 
diagrammatic axiomatization and then we study the Hilbert representations, where graph- 
theoretic objects shall be associated with concrete elements of Hilbert spaces. 

3.1. Diagrams and tensors. 

Diagrams, boxes, decorations and wires. Our starting point is a set S endowed with 
an involution without fixed point *. The set S splits as S S* according to the involution. 
Elements of S are called decorations. 

A diagram is a collection of decorated boxes and possibly wires (or strings) connecting the 
boxes along their decorations according to rules which we shall specify. In terms of graph 
theory, a diagram is an unoriented (multi-)graph whose vertices are boxes, and whose edges 
are strings. Each vertex comes with a (possibly empty) n-tuple of indices (or decorations or 
labels) in S". The number n of decorations may depend on the vertex. We say that two 
diagrams are isomorphic if they are isomorphic as multi-graphs with labeled vertices. 

A box is an elementary diagram from which we can construct more elaborate diagrams 
by putting boxes together and possibly wiring them together. Each box i? of a diagram has 
attached to it a collection of n{B) decorations in S'^^^\ The union of the decorations attached 
to a box B is denoted by S{B) U S*{B). 

Graphically, boxes are represented by rectangles with symbols corresponding to the deco- 
rations attached to them (see Figured]). We take the convention that decorations in S* are 
represented by empty (or white) symbols and decorations in S by full (or black) symbols. 
Each decoration is thought as having potentially up to two attachment points. An inner one 
(which is attached to the box it belongs to) and an outer one, which we shall allow to be 
attached to a string later on. 



I 

I M » 



Figure 1. A box M 
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Constructing new diagrams out of old ones. Given a family of existing diagrams (e.g. 
boxes) there exists several ways of creating new diagrams. 

(1) One can put diagrams together, i.e. take their disjoint union (when it comes to 
taking representations in Hilbert spaces, this operation will amount to tensoring). 
One diagram can be viewed as a box. This amounts to specifying an order between 
the boxes. 

(2) Given a diagram A and a complex number x, one can create a new diagram A' = xA. 

(3) Given two boxes A, B having the same n-tuples of decorations, one can define A + B. 
This axiom and the previous one (together with evident relations such ss, A + A = 2 A 
which we don't enumerate in detail) endow the set of identically decorated diagrams 
with a structure of a complex vector space. 

(4) One can add wires to an existing diagram (or between two diagrams that have been 
put together). A wire is allowed between the outer attachment of two decorations 
only if the decorations have the same shape and different shadings. Such a wire can 
be created if and only if the two candidate decorations have their outer attachments 
unoccupied. 

(5) There exists an anti-linear involution on the diagrams, denoted by *. This operation 
does nothing on the wires. On the boxes, it reverts the shading of the decorations. 
The involution * is conjugate linear. 

Hilbert structure. We shall now consider a concrete representation of the diagrams intro- 
duced above as tensors in Hilbert spaces. We start by assuming that the set S of full (or black) 
decorations corresponds to a collection of finite dimensional Hilbert spaces 5* = {Vi, V2, • • • }. 
An important fact that will be useful later is that each Hilbert space Vi comes equipped with 
an orthonormal basis {ei, 62, ... , Cdim Vil- Our aim is to define a *-linear map T between the 
diagrams and tensors in products of Hilbert spaces in the above class and their duals. By 
duality, white decorations correspond to dual spaces S* = {^1*, V'2*, • • • }. With these conven- 
tions, boxes can be seen as tensors whose legs belong to the vector spaces corresponding to 
its decorations. In a diagram, symbols of the same shape denote isomorphic spaces, but the 
converse may be false. A particular space Vi (or V*) can appear several times in a box. The 
reader acquainted with quantum mechanics might think of white shapes as corresponding to 
'bras' and black shapes corresponding to 'kets', but we shall get back to quantum mechanical 
notions later. 

To a box B we therefore associate a tensor 



(7) 


Tb g 


(g) V 




(g) v; 






_i€S{B) 




jes-iB) 



Using the canonical duality between tensors and linear applications, Tb can also be seen as a 
linear map 

Tb: (g) Vj^ Vi, 

jes*{B) ieS{B) 

We use freely partial duality results, and for example, an element of V <^ W* can as well be 
seen as an element of C{W, V) or C{V*, W*). 

Equation ([7]) defines the map T from the collection of boxes to the collection of vectors in 
Hilbert spaces obtained by tensoring finitely many copies oiVi,iG S{B) U S*{B). This map 
is denoted by 

T : B^Tb 
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and we now explain how we can extend it to all diagrams. A wire connecting two decorations 
of the same shape (corresponding to some Hilbert space V) is associated with the identity map 
(or tensor) 1 : V ^ V. Together with our duality axiom, it also corresponds to a canonical 
tensor contraction (or trace) 

C -.V* (S)V ^C. 

We denote the set of wires in a diagram V by C{p). 

With this notation, a diagram T> is associated with the tensor T obtained by applying all 
the contractions ("wires") to the product of tensors represented by the boxes. One is left with 
a tensor 



Tb 



\B box of V 



n c 

cec{v) 

This is well defined (provided that one specifies one total order on the boxes): the order of 
the factors in the product does not matter, since wires act on different spaces. For a box B, 
we denote by FS{B) C S{B) the subset of black decorations which have no wires attached 
(we call such a decoration free). FS*{B) is defined in the same manner for white decorations 
(dual spaces). With this notation, the tensor T-d associated to a diagram V can be seen in 
two ways: as an element of a Hilbert space 



V* 







Vi 
Vi. 



or, equivalently, as a linear map 

We need two further axioms to ensure that we are indeed dealing with acceptable Hilbert 
representations. 

(1) A diagram such that all outer attachments of its decorations are occupied by wires 
corresponds canonically to an element in C. In addition, a trivial box with a given 
decoration of type i closed on itself by a wire into a loop takes a value in N. This 
value is called the dimension of Vi. 

(2) Given a diagram P, if it is canonically paired to its dual V* by strings, the result lies 
in M+. 

Special diagrams. To make our calculus useful, we need to introduce a few special diagrams 
(equivalently, boxes) satisfying some specific axioms. 

(1) The trivial box. A wire connecting two identically shaped decorations of different 

shading corresponds to the identity map \:V ^V. We shall call this box the trivial 

or the identity box. 



Figure 2. Trivial box 
It satisfies the following identity axiom: 



Figure 3. Trivial axiom 
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(2) Bras and kets. The simplest boxes one can consider are vectors and linear forms. 
Following the quantum mechanics 'bra' and 'kef vocabulary, vectors, or (1, 0)-tensors 
have no white decorations and only one black decoration, whereas linear forms (or 
(O,l)-tensors) have one white label and no black labels. Since our Hilbert spaces come 
equipped with fixed basis, we introduce some special notation for the ket |ei) = ei 
and the bra (ei| = (ei, •) corresponding to the first vector, as in Figured] (the choice 
of the first element of the basis being of course arbitrary) . 



o = c ei 



Figure 4. Fixed ket and bra 

(3) The Bell state. Since each space V £ S comes equipped with a particular fixed basis 
{ej}^l™^, we can define the bra Bell state as the tensor (it is in fact a linear form) 

dim V 



BelVy = 



and its ket counterpart (which is a vector mV<SiV) 

dim V 

Belly = ei® ei. 

i=l 

This notation is needed in the sense that Bell states are not canonical and are not 
well defined from the sole data of V. They rely on some additional real structure 
of the vector space V which can be encoded by the data of an explicit basis. Bell 
states are represented in Figure 5(a) They satisfy the graphical axiom in Figure [5(b)J 



Bell states play a central role in our formalism; we shall see later that they allow us 
to define the transposition of a box and even to consider wires connecting identical 
decorations. 




(a) (b) 
Figure 5. Bell states and axiom 



(4) Unitary boxes. Boxes associated to unitary matrices U satisfy the graphical axiom 
depicted in Figure [6] which corresponds to the identities UU* = U*U = 1. 



c U 



U* 



c u* 



u 



Figure 6. Unitary axioms 
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M 
(a) 



M p-^ y 

(b) 

Figure 7. Some simple diagrams 



N 
(c) 



3.2. Examples. Let us now look at some simple diagrams which illustrate this formalism. 

Suppose that each diagram in Figure [7] comes equipped with two vector spaces Vi and V2 
which we shall represent respectively by circle and square shaped symbols. In the first diagram, 
M is a tensor (or a matrix, depending on which point of view we adopt) M £ V{ <^ Vi, and 
the wire applies the contraction f^iVi C to M. The result of the diagram T>a is thus 
Tp^ = Tr(M) G C. In the second diagram, again there are no free decorations, hence the 
result is the complex number Tx)^ = {y,Mx). Finally, in the third example, is a (2,2) 
tensor or a linear application G L(Vi V2, Vi V2). When one applies to the tensor 
the contraction of the couple {Vi, V^), the result is the partial trace of N over the space Vi: 

Bell states allow us to introduce the transposition operation for a tensor (or a box) as 
follows. We define transposition for a matrix M (or a tensor M € V* V) and we extend 
it in a trivial way to more general situations. Graphically, the box corresponding to the 
transposed tensor M* is defined in Figure 8(a) it consists in connecting an appropriate Bell 



state to each decoration of the box. Note however that this operation is different from the 
involution * applied to the same box. Moreover, Bell states allow for wires connecting identical 
shaped symbols of the same color, as in Figure [8(b)[ Such non-canonical tensor contractions 
{V (Si V ^ C or V* V* ^ C) are shorthand graphical notations for the corresponding 
diagram containing a Bell state, and we shall use them quite often in what follows. 



c M* • = 



(a) 




Figure 8. Bell states and transposition 



Also, for reasons which shall be clear later, we shall sometimes make substitution M = 
(Af*)*. Finally, by grouping two Bell states together, one obtains the (non-canonical) tensor 
E (Figure [9]), called "the maximal entangled state". It corresponds to the tensor 

dim V dim V 

E= J]] e* e* (g) Cj ^Cj £V* 0V* ^V^V. 
i=i j=i 

The reader with background in quantum information will notice that the maximally entangled 
state we just defined is not normalized in order to be a density matrix. The reader with 
background in planar algebra theory will recognize a multiple of the Jones projection. 

The diagram in Figure [10] is of crucial importance in what follows. It corresponds to a 
quantum channel in its Stinespring representation (see Eq. ([6])). Round shaped inputs and 
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Figure 9. The maximal entangled state E 



outputs correspond to the Hilbert space TC and squares correspond to the "environment" 
We shall study such diagrams with random interaction unitaries U in Sections [5] and El 



/C. 



c • 



U 



Y 



U* 



Figure 10. Diagram for a quantum channel in its Stinespring form 



3.3. Comments on other existing graphical calculi. The above formalism is the one 
that seemed the most compatible with Weingarten calculus. Here, we comment about already 
existing graphical formalisms, in the hope that this section will serve as a dictionary for the 
reader acquainted to one of the calculi below. 

Our calculus is mainly inspired by Bob Coecke's Kindergarten Quantum Mechanics [1]. 
However we choose not to orient the strings; rather, we separate with color (black/white) the 
vector spaces and their duals, therefore there is only one possible pairing. A common feature 
of the two calculi is the central place occupied in the formalisms by Bell states. 

V.F.R. Jones's theory of planar algebras [10] is also connected to our graphical calculus. 
One of our diagrammatic axioms is the existence of a Bell state. This is very closely related 
to the axioms of Temperley Lieb algebras and the diagrammatic for a Jones projection. Most 
of our calculus could take place in Jones' bipartite graph planar algebra. 

4. Planar expansion 

In this section, we consider diagrams that may contain random matrices. This is where 
probability theory comes into play; we focus on the case where the random elements appearing 
in the diagrams are random unitary Haar-distributed on some finite dimensional unitary 
group. 

Our task is to compute expectation values E(D) of diagrams T) containing boxes associated 
with random unitary operations. We shall write E(2?) as a weighted sum of some diagrams 
obtained from T) that to not contain anymore random tensors. 

4.1. Expectation of a diagram containing random independent unitary matrices. 

Suppose we have a diagram V that has boxes of two types: either boxes of type U, U*, U or 
[/* where C/ is a unitary random variables in a fixed space of type End((8)jg/Vi), distributed 
according to the Haar measure on the unitary group U{n) on this space, or other boxes which 
are independent (as classical random variables) from U (this includes deterministic boxes or 
tensors). We shall now present an algorithm for computing the expectation of such a diagram 
with respect to the probability law of U. Before describing the algorithm, let us note that if 
a diagram contains several independent Haar unitary matrices, one can recursively apply the 
algorithm to compute the expectation over all the random unitary matrices appearing in the 
diagram. 
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The first step in our algorithm is to ensure that D contains only boxes of type U and U. 
This can be done by using the transposition transformation via Bell states, and replacing [/* 
boxes by U boxes with the opposite shading of the decorations and U* boxes by U boxes. 

Next, we introduce a concept of removal of boxes U and U. A removal r is a procedure 
which transforms a diagram D into a new diagram Dr which does not contain neither U nor 
U boxes. In other words, r is a way to remove random unitaries U from a diagram D. The 
set of all admissible removal procedures for a diagram D will be denoted by Rem(!D). 

We now move on to describe removal procedures and how they operate on diagrams. First 
of all, a removal is not possible if the number of boxes U in D is different from the number 
of boxes U. In such a case, the set Rem(I?) will be defined to be empty. This rule is the 
diagrammatic equivalent of Eq. ([2|) from Theorem 12.21 

Assuming that the number of U boxes and U boxes is the same, a removal r is a way to 
pair decorations of the U and U boxes appearing in a diagram. More precisely, r is the data 
of a pairing a of the white decorations of U boxes with the white decorations of U boxes, 
together with a pairing /3 between the black decorations of U boxes and the black decorations 
of U boxes. Assuming that V contains p boxes of type U and that the boxes U (resp. U) are 
labeled from 1 to p, then r = {a, (5) where a, (5 are permutations of Sp. 

Given a removal r G Rem(!D), we construct a new diagram Dr associated to r, which 
has the important property that it no longer contains boxes of type U or U. We proceed 
in the following way: one starts by erasing the boxes U and U but keeps the decorations 
attached to them. Assuming that one has labeled the erased boxes U and U with integers 
from {1, . . . ,p}, one connects all the (inner parts of the) white decorations of the i-th erased 
U box with the corresponding (inner parts of the) white decorations of the a(z)-th erased 
U box. In a similar manner, one uses the permutation /3 to connect black decorations. A 
diagrammatic explanation of the above algorithm is described in Figure [TTl 



I ~ 
^ [7 Ha 



a{i) ^ [/ P 

Figure 11. Removal: elimination of boxes and pairing of decorations 

We are now ready to state the main result of this section. 
Theorem 4.1. The following holds true: 

Eu(V)= Yl VrWg{n,ar^). 

r={o,/3)GRem(X>) 

Proof. This is just Weingarten calculus of Theorem 12.21 applied to our graphical conventions. 

□ 

When more than one independent unitary matrices U,V, . . . are present in a diagram, we 
proceed by induction over each independent matrix: we successively remove U, V, etc. One 
can check directly that the order of the induction does not change the final result. This is 
compatible with the probabilistic property of the expectation, E[/y(D) = Ev/(E[/(D)). 
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Figure 12. The diagram for Tr(yP) 



Theorem 14.11 might just look as a reformulation of known results, but without this graph- 
ical method, obtaining the main results of this paper is extremely cumbersome and very 
counterintuitive. 

Let us now comment on the first step of our removal algorithm, replacing U* boxes with U 
boxes. The purpose of such a substitution is purely practical: later in the removal procedure, 
we pair decorations of the same color. If we should have decided to work with U and U* 
boxes, one should always pair decorations of different colors, and this can turn out to be 
rather cumbersome when doing combinatorics. On the other hand, each time we replace a U* 
box by a ?7 box, we introduce two more Bell states into our diagram (see Figure [8|); although 
we decided not to display such states and rather to allow wires connecting decorations of the 
same color, this operation increases in some sense the "complexity" of the diagram. 

In the next sub-section we present a warm-up toy example of Theorem 14.11 Further appli- 
cations of the above theorem will be considered in a forthcoming paper where problems 
from free probability theory will be treated using similar techniques. 

4.2. First example: partial tracing a randomly rotated matrix. As a first application 
of the graphical formalism, we consider the following problem. Let X € A4„jt(C) be a de- 
terministic matrix. In a manner similar to random quantum channels, we define, for a fixed 
integer parameter k ^ 1, the random matrix 

Y = Tvk[UXU*]£ Mn{C), 

where U G U{nh) is a Haar distributed random unitary matrix. Notice that we are considering 
here non-normalized traces. Using our graphical formalism, we shall compute the moments of 
y, E[Tr(y^)] for all p ^ 1. After replacing U* boxes with U, one gets the diagram in Figure 
WH where round decorations correspond to C" and square ones to C'^. Note that in Figure \T2\ 
there are p groups ^^UXU" wired together. 

By Theorem 14. 1^ the expectation (with respect to the Haar measure of the unitary group 
U{nk)) of this diagram is a weighted sum (with Weingarten weights) of diagrams T>r obtained 
after the removal of U and U boxes. Such diagrams Dr contain only X blocks and loops of 
different types. Let us compute the the value of a diagram P^, where the removal is given by 

Suppose we number the boxes from 1 to p and the permutations a connects the white 
decorations and P the black decorations. The point here is that the permutation a will be 
responsible for the loops appearing in T>r, whereas P will connect X boxes. We start by 
counting the number of loops in Dr- Loops are of two types, the ones containing the top, 
square decorations (a loop of this type has a value of k) and the ones that come from the 
bottom, round decorations, each with a value of n. Since the top decorations are already 
connected (in the original diagram D) by the identity permutation, the number of loops of 
the first type is given by the number of cycles of a, they give a total contribution of A;*°. 
Loops corresponding to round decorations are initially connected by the cycle 



(8) 



7 = (n n - 1 • • • 3 2 1) G 5p. 
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Hence, the number of such loops is #(7 ^a) and they give a total contribution of n*^''' In 
conclusion, the total contribution of loops is k'^'^n'^^"' "\ The contribution of the X boxes 
is straightforward to compute, since these boxes are connected only by /?. After the removal 
we get #/? connected components of powers of X, and the total contribution is Tr^(X). 
Putting all this factors together, we obtain the following proposition: 

Proposition 4.2. The mean p-th moment of the random matrix Y = Tr^.[?7XC/*] is given by 

(9) E[Tr(yf)]= A:#"n#(^"'")Tr^(X)Wg(a/3-i). 

In the particular case where X is a rank one projector, one has Tr/3(X) = 1 for all permu- 
tations /3 and since 

/p-i 

J2 Wgink,a) = Ylink+j) 

a£Sp \j=0 

one obtains the following simplification: 

(10) E[Tr(yP)] = I J|(n/t + i) I ^ A:#"n#(^"'"). 

\i=o / aeSp 

We will discuss at length the simplifications that occur when dealing with rank one projectors 
in the forthcoming paper [4]. 

5. Tensor products of independent random quantum channels 

In the present section and in the next one, we consider two different models of tensor 
products of random quantum channels. In both cases, we first fix an interesting input state 
X12 (the Bell state) and investigate the random matrix 

[$10«>2](Xi2). 

Our two models correspond to the choice of two different random pairs ($i,<I>2)- In both 
models, the channels ^>i, ^2 are defined by random unitary matrices Ui, U2 via the Stinespring 
representation introduced in Eq. Q. 

The difference between the two cases lies in the correlation between the random matrices 
Ui and U2- In the first model, interaction unitaries are independent Haar distributed 
random matrices. The second model, which is more involved, deals with the case where Ui is 
distributed according to the Haar measure and U2 = Ui. This choice introduces an interesting 
symmetry into the problem and such a construction has become classical in quantum infor- 
mation theory [TJ |6] . Asymptotic results in this case shed light on the interesting phenomenon 
that the output of the product channel has one "large" eigenvalue. We call this phenomenon 
the Bell state phenomenon. In order to simplify the notation, we shall assume that Yi and Y2 
are rank-one projectors and that ni = n2 = n, fei = /c2 = k. 

Before looking in detail at the two models of interest, let us make one brief comment on 
the choice of the input state of the channels. It is clear that if one chooses an input state 
which factorizes X12 = Xi X2 , then 

[$1 (g> ^2]{Xi2) = ^l{Xi) (g> $2(^2), 

and there is no correlation (classical or quantum) between the channels. In order to avoid 
such trivial situations, one has to choose an input state which is entangled. An obvious choice 
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(given that ni = n2 = n) is to take X12 = En, the n-dimensional Bell state (or the maximally 
entangled state), and we shall use this state in what follows. 

The first model, although new, does not bring strikingly new results from the random 
matrix point of view. We treat it here as an illustration of what our calculus can allow to 
compute, and as a point of comparison with the second model. 

Independent interaction unitaries. In the remaining of this section, we consider two 
independent realizations Ui = U and U2 = V oi Haar-distributed unitary random matrices 
on U(nk). For both channels the state of the environment is a rank-one projector and we are 
interested in the x random matrix 

Z=[cI>^®<I.^](i?„), 

where En is the maximal entangled Bell state (notice the 1/n normalization) 



En 



1 " 



\ei){ej\. 



The diagram associated with the (2,2) tensor Z is drawn in the Figure [131 we chose to 
represent by squares decorations corresponding to C'^ and by circles decorations corresponding 
to C". 



z = $^ ®$^(i;„) 



c 

U 

c 


■— □ 1 ■ 

n 

• — Ck J» — 


I 

U* 

» 








c 

V 

c 


• — Of N» — 
■— □ ■ 


V* 

I 



Figure 13. Z = (g) $^(^„) 



As usual, we are interested in computing the moments E[Tr(ZP)] for all p ^ 1 using the 
graphical method. We start by replacing U* (resp. V*) blocks by U (resp. V) blocks. An 
important point here is that there are two type of blocks corresponding to the independent 
random unitary matrices U and V (when computing the p-th. moment of Z, there are p blocks 
of each type). This has two important consequences: when expanding the diagram in order 
to compute the expectation of the trace, one can only pair U blocks with U blocks and V 
blocks with V blocks; "cross-pairings" between U blocks and V blocks are not allowed by 
the expansion algorithm. This algorithm proceeds iteratively, first by removing, say, the U 
blocks (the V blocks being treated as constants) and then by removing the V blocks. Hence 
"cross-pairings" cannot occur. The second consequence of the presence of two independent 
Haar unitary matrices is that in the final expression for the expectation of the diagram, there 
will be two Weingarten weights, one for each independent unitary integration. 

Lemma 5.1. The following holds true (7 is the cycle permutation defined in Eq. 

(11) E[Tt{ZP)]= A;#"t'+#'^^n#(^"'°^)+#(^"'°^)+#(^^''^^)~PWg(ac7/3j;i)Wg(av'/3^i). 
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Proof. As it has already been stated, the expectation with respect to both unitary matrices 
U and V can be seen as the result of two removal procedures, and hence the Weingarten sum 
shall be indexed by a pair of removals {ru,rv)- In other words, the Weingarten sum shall 
be indexed by 2 pairs of permutations, one for each type of block; we shall denote them by 
o^u, /5c/, ay, Pv £ 'Sp- The four permutations are responsible for pairing blocks in the following 
way (1 ^ i ^ p): 

(1) the white decorations of the i-th ?7-block are paired with the white decorations of the 
au{i)-th U block; 

(2) the black decorations of the i-th U -block are paired with the black decorations of the 
Pu{i)-th. U block; 

(3) the white decorations of the i-th. F-block are paired with the white decorations of the 

Q;v'(i)-th V block; 

(4) the black decorations of the i-th y-block are paired with the black decorations of the 
/3y(i)-th V block. 

The diagram associated with Tr(Z^') contains, aside from the random unitary blocks, deter- 
ministic bras □ and kets ■. However, these boxes have a trivial contribution of 1 to the final 
result. Hence, the result of the graph expansion is a (sum over a) collection of loops, multi- 
plied by some scalar factor. The different contributions of a quadruple (»(/, ay, G 
are given by (recall that circles correspond to ra-dimensional spaces and squares correspond 
to fc-dimensional spaces): 



(1) 


"□i7" -loops: 




(2) 


"□F" -loops: 




(3) 


"o?/" -loops: 


^#(7"^o;[/). 


(4) 


"oy"-loops: 




(5) 


"■t/"-loops: 


none; 


(6) 


"■F" -loops: 


none; 


(7) 


"•?7,y" -loops: n*^^u M 



(8) normalization factors 1/n from the Bell matrices En- n~^; 

(9) Weingarten weights for the i7-matrices: Wg{auP^^)■, 
(10) Weingarten weights for the F-matrices: Wg{av/3y^). 

Adding all these contributions, we obtain the announced exact closed-form expression. □ 

Asymptotics. The preceding expression is intractable at fixed n and k, so we study two 
asymptotic regimes: 

(I) n fixed. A; — > oo; 
(n) k fixed, n — > oo; 

At this stage, before looking into each particular asymptotic regime, we can make an im- 
portant observation. Notice that in the preceding expression, aside from the factor rC^'^^u ^^), 
the general term of the sum factorizes into a "(«[/, /3t/)" part and a "(ay, /3y)" part. This was 
to be expected, since the coupling between the two channels is realized by the input state 
which has a contribution of rt^^^u Pv)-p _ Let us also note that a third interesting asymptotic 
regime, n, A;— >oo, A;/n— >c>0 will be studied using the same methods in a forthcoming 
paper. 
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Theorem 5.2. In the firs regime, n fixed, k — > cxd, the output of the tensor product of the 
channels, for large values of k, is close to the chaotic state 

p* = —■ 

In the second regime, k fixed, n — > oo, the asymptotic eigenvalues of Z are 1/k'^ with multi- 
plicity k"^ and with multiplicity — k"^ . 

Proof. Using the standard asymptotics for the Weingarten functions 

WgiauPu^) {nk)-'P'\'^uPd^\ Moh[auPu^) and 
Wg(ay/3^i) {nk)-P-\"^^v\ Moh{avl3y^), 
the power of k appearing in a general [au , f3u , ay , j3v) term is 

jj)-\au\+p-\av\-p-\auP^^\-p-\avPv^\ — J^-{Wu\ + Wv\ + \o^U Pu^\ + \<^v Pv^\) 

It is obvious that all the terms converge to zero, except the one with au = Pu = otv = Pv = id. 
Using Mob (id) = 1, we conclude that 

lim E[Tr(ZP)] = n^'^P. 

One can restate this in terms of the empirical eigenvalue distribution of the x matrix 
Z: 

1 

i=l 

In other words, the output of the tensor product of the channels, for large values of k, is close 
to the chaotic state 

P* = —■ 

As for the second regime, using similar considerations, we obtain 

lim E[Tr(ZP)] = fe^-^P. 

n 

Note that both regimes presented here are trivial to some extent. We could prove at only 
a small additional cost that the convergence of the eigenvalues is actually almost sure. See 
the next section for the technique of proof, a direct adaptation of the Borel-Cantelli lemma. 

Finally, we want to emphasize that the asymptotic behavior of the output in the second 
regime changes drastically when considering the conjugate case, and that this fact will have 
very important consequences in the theory of quantum information. 

6. Tensor products of conjugate random quantum channels 

We have seen that tensor products of independent random channels have an eigenvalue 
behavior close to the single channel case (see [Hj for the treatment of the single channel case) 
- despite the fact that the input state is maximally entangled. In this section, we consider 
the case where Ui = U , U2 = U and U G lA{nk) is a Haar uniform random unitary matrix. 
Tensor products of channels of this type are now classical in the literature (see [3|6]). One of 
the reasons why channels of this particular form receive such attention is that one can show 
that the product channel has a "trivial large eigenvalue" of order 1/k. We shall provide a 
graphical proof of this fact later, in Lemma 16.61 
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Again, we are interested in the moments of the random matrix Z = ^ 0^ (En), depicted 
in the Figure [Ml 



u 



u* 



» — Ck /• — c 

• — cr ^»—c 



U 



U 



u* 



u 



u 



u 



3 



Figure 14. Z = ® 



This time calculations are more complicated, because only one unitary matrix appears in 
the product channel. This means that in the removal algorithm, one can pair boxes from $^ 
with boxes from thus obtaining more complicated patterns. Another consequence of the 
fact that we use only one random unitary matrix is that the Weingarten sums are indexed by 
only one pair of permutations (a, f3) £ ^fp- 

In order to count the loops obtained after the graph expansion, we label the U and the U 
boxes in the following manner: 1-^, 2^ , ■ ■ ■ for the U boxes of the first channel (T as "top") 
and 1^, 2-^, . . . for the U boxes of the second channel (B as "bottom"). We shall also order 
the labels as {l'^, 2^, . . . , 1^,2^,... ,p^} ~ {1, . . . , 2p}. A removal r = (a, /3) G 5|p of the 
random (U and U) boxes connects the decorations in the following way: 

(1) the white decorations of the i-th [/-block are paired with the white decorations of the 
a(«)-th U block; 

(2) the black decorations of the i-th [/-block are paired with the black decorations of the 
/?(z)-th U block. 

Next, we introduce two fixed permutations 7, (5 G S2p which will be useful in counting the 
loops. The permutation 7 represents the initial wiring of the o decorations (before the graph 
expansion) and 6 accounts for the wires between the • decorations (which come from En). 
More precisely, for all i, 

(12) ^{f) = {i-lf, = and 6{f) = i^, 6{i^)=f. 
We are now ready to compute the average moments of the random matrix Z. 

Lemma 6.1. The following holds true 

(13) E[Tt{ZP)]= A:#"n#("^'')+#(^^)-PWg(a/3-i). 

Proof. With the notations introduced above, we can now count the contributions for each 
individual pairing (a, /3): 

(1) "□"-loops: 

(2) "o"-loops: n#('^'^''); 

(3) """-loops: none; 

(4) "•"-loops: n^^^l^^ ) = n^(f^^) (notice that 6 is an involution); 

(5) normalization factors 1/n from the En matrices: n~P; 

(6) Weingarten weights for the [/-matrices: Wg{af3~^). 

Adding up all the above contributions, we obtain the claimed formula. □ 
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In the rest of the section, we shall focus on the asymptotic regime n fixed, k ^ oo, and in 
the next section we shall look into the more interesting case k fixed, n ^ oo. 

Before stating the asymptotic result, let us make two preliminary remarks. In the case of 
the conjugate product channel, since only one unitary matrix appears in the diagrams, there 
is a notable difference concerning the asymptotics of the Weingarten function: 

Wg{ap-^) ~ (nA:)-2p-l°/3~'l Mob(a/3~^). 

One can easily compute the following quantities involving the permutations 7 and 6 which 
will be useful later, when doing asymptotics: I7I = 2p — 2, \6\ = p, |7(5| = p. 

Proposition 6.2. In the asymptotic regime where n is fixed and — > 00, the random matrix 
Z converges to the chaotic state 

P* = 

Proof. Computing the asymptotic trace for large k gives 

E[IV(ZP)]~ Yl A:-(l"l+l"^"'l)nP-(l"^"'l+l^^l+l"'3"'l)Mob(ar')- 

Minimizing the power of k above gives \a\ + |a/3~^| ^ 0, with equality iff a = /? = id, hence 

lim E[Tr(ZP)] = n^-^P, 

and the conclusion is the same as in the case of two independent quantum channels: the 
output Z is asymptotically close to the chaotic state □ 

6.1. Conjugate channels, the Bell phenomenon. We are left with studying the most 
interesting regime, k fixed and n 00. Our main result is as follows: 

Theorem 6.3. In the regime of k fixed, n 00, the eigenvalues of the matrix Z converge 
almost surely towards: 

• ^ + — p-, with multiplicity one; 

• p- — with multiplicity k"^ — 1; 

• 0, with multiplicity — k"^ . 

Note that it follows from the Stinespring theorem that Z has at most k'^ non-zero eigen- 
values. Therefore a moment approach is possible for the proof. 

We start with a technical Lemma about the structure of geodesies between the specific 
permutations 7 and 5 introduced in Equation ()12p . 

Lemma 6.4. For 1 ^ i ^ p, let Ti he the transposition {{^ , [i — 1)'^) . Then the permutations 
a on the geodesic 7 — > a — > 5 are indexed by subsets A CI {1, . . . ,p}: = 7 Higa '^i- Moreover, 
for such a permutation, we have 

\2p-2 if A = 9, 
[2p-\A\ if A ^9. 

Proof. If ^ = 0, then \a\ = \ j\ = 2p — 2. Otherwise, after computing the action of a 



\a\ 



rp \i^ if i G A, 



a{i^) 



(i-lf Ifi^A; 

'i^ if (i + 1) G A, 

(i + 1)^ if (i + 1) ^ A; 
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it is easy to see that each element i A spans a cycle of a and thus |a| = 2p — = 
2p-\A\. □ 

We split this proof of Theorem 16.31 in two steps: first we prove the convergence in expecta- 
tion, and then we prove the almost sure convergence. 

Proof of the convergence in expectation. Using the same asymptotic formula as in the previ- 
ous section (this time for large n), the quantity one wants to minimize in this case is 

|a7-^| + \P6\ + \a(3-^\ = \-f-^a\ + la"^/?] + \(3~^d\ ^ \-f-^6\ = p. 

Equality is attained when 7^a^/3^(5isa geodesic in S2p- Using this observation, we 
obtain 

E[T¥(ZP)]~ A:-(l"l+l°^''l)Mob(a/3-^). 

It turns out that we can compute exactly the last sum as follows. First, notice that the 
geodesic condition 7 ^ q 5 can be restated as id — > j~^a — > 7~^/? — > j~^5. But 

7~^5 is a product of p transpositions with disjoint support: j~^5 = HiLi'^i' where is the 
transposition (z^ , (i — 1)^) for all 1 ^ i ^ p. As in Lemma 16.41 permutations on a geodesic 
between id and j^^5 are parameterized by subsets of {1, . . . ,p} as follows. Permutations 7~^a 
and 7~"^/3 lie on a geodesic between id and 'y~^d (i.e. id 'y~^a — > 7~^/? — > 7~^5) if and only 
if there exist two subsets C A C B C {1,2, ... ,p} such that 

For two such permutations, it is obvious that |a~^/?| = |(7^^q;)^^7~^/3| = |i? \ A|. In order 
to compute |q:|, we rely on Lemma |6.4[ 

Since a~^P is a product of |i? \ A| transpositions of disjoint support, it follows by ^5] that 
Mob(Q-V) = (-l)l^\^l and we are left with the following expression: 

E[TV(Zf)]~^A;-(2p-2+|B|)(_i)|s|+ ^ k-'-^P-\A\+\B\A\)^_-^^\B\A\_ 

B ID^ACB 

Using the multinomial identities 

x'^l = (1 + xf and 
0cAc{i,...,p} 

Y = (l + x + yf, 

0CACBC{l,...,p} 

we obtain the asymptotic traces of the output matrix Z: 

We conclude that the matrix Z has, asymptotically, the following eigenvalues: 

• 1 + - -g^, with multiphcity one; 

• ~ ■p! with multiplicity k"^ — 1; 

• 0, with multiplicity n'^ — k'^. 

□ 
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Next we move on to the proof of almost sure convergence. We would like to mention about 
the proof below that we believe that it should be possible to prove that 

E Utt{ZP) -ETT{ZP)f] = 0(n-2) 

simply by observing that the function 

U Tt{ZP) 

is Lipschitz on the unitary group and by applying a Gromov-Milman type concentration mea- 
sure argument. We refer to [8j for an exposition of such techniques. The authors acknowledge 
that this approach might be slightly less cumbersome in the specific case of this proof. How- 
ever, we chose to keep our proof of a combinatorial nature for the sake of coherence. 

Proof of the almost sure convergence. It is a standard technique in probability theory that in 
order to show the almost-sure convergence of the eigenvalues of Z to their respective limits, 
it suffices for the covariance series to converge, for all values of p: 



[(Tr(ZP) -ETr(ZP))2 

n=l 



< OO. 



Indeed, this inequality together with the Borel-Cantelli lemma imply that almost surely as 
n — > oo, 

Tr(ZP) ^ETr(ZP) 

The two ingredients which make the proof work are the following. The first one is the fact 
that the error one makes when approximating the Weingarten function with its dominating 
term is of the order —2: 

Wg(a) = (nA:)-(P+l°l)(Mob(Q) +0((nA:)-2)). 

This follows from Theorem 12.31 and the definition of Mob below. The second ingredient is 
contained in the geodesic inequality |7~"^a| -|- |a~^/?| -|- l/?""*^^] ^ |7~"^<^| = P- Earlier, we have 
completely described the set of couples (a, /3) which saturate the equality. It turns out that 
one can say more on the values of the function (a, /?) i— > E{a, j3) = |7~^a| -|-|a~^/3| -|-|/3~^(5| —p, 
as follows. Applying Lemma 12.41 two times, it is clear that the values taken by E{a,(5) are all 
even, E{a., 13) ^ 1 and thus 



E[Tr(ZP)] = A;#°n#("^ ')+#(^^)-PWg(a/?-^) 



E 



^-(|a|+|a/3-i|) Mob(a/3-i) + o(\ 



We start by computing the easiest term of the covariance, namely 
E[Tr(ZP)]2 = Y fc"^'°''+'°^^i~''+'"''+'°'^2'''^Mob(ai/3fi)Mob(a2/32~^) + . 

The second term E[Tr(ZP)^] is more difficult to estimate and one needs to introduce the 
permutations 7, 5 S : 

7 = (1^2^ • ■■p^){{p + Ifip + 2f--- {2pf){l^2^ ■ ■■p^){{p + Ifip + 2f... {2p)% 
-5 = (1^1^)(2^2^) . . . (pV)((p + lYiP + If) ■ ■ ■ {{2p)^{2p)^). 
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With this notation, we have 

^ ^-(|„| + |„/3-i|)^2p-(|7-la| + |a-l/3| + l/3-l<5|) Mob(a/3-') + O 



1 



One can easily show that |7~^a| + + |/3 ""^^l ^ It"^*^! = 2p. Since both 7 and 5 



leave invariant the sets {1 



T,B r,T,B 



. . . ,p^'^} and {{p+lY'^, + . . . ,2p^'^}, geodesic 



couples (a, /5) are obtained as direct sums 

a = ai © 02, 

where ai,/3i G 2^'^, . . . a2p2 G + if (p + 2f . . . , 2p'^'^}) are 

such that 7i — > CKi ^ /?! — > 5i and 72 ^ 02 ^ /32 — ^ ^2 are geodesies (the permutations 71^2 
and 61^2 are defined in an obvious way). One has also that \a\ = \ai\ + \a2\, \P\ = + \P2\ 
and that Mob(a/3~^) = Mob(ai/?f^) Mob(a2/3^^). Putting all this together, one gets the final 
expression which ends the proof 

E[TV(ZP)2] = ^-(|ai| + |ai/3-i| + |a2| + |a2/3-M) Mob(ai/3f 1) Mob(a2/?2"^) + ^ ( ^) ■ 



n 



6.2. Generalization of Theorem 16.31 and discussion. We finish this paper by studying a 
generalization of the model investigated in the previous section. We consider the case where 
k is a fixed integer, and t £ (0, 1) is a fixed number (possibly a function of k). For each n, we 
consider a random unitary matrix U G A4„fc(C), and a projection qn of A4„fc(C) of rank pn 
such that pn/{nk) ~ t as n — > cxd. Our model of a random quantum channel is 

<^ : Mp„{C) ^ Mn{C) 

given by 

$(X) = Tr,,(;7X;7*) 

where the density matrix X satisfies X ^ qn (in other words we consider the isomorphism 
Qn-M.nk{'C)qn — A^p„(C)). Graphically, our model amounts to Figure [T5l 



> $^(£;t„fe) 





k-A ^_ A-t. 






tnk 


u* 
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U* 
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u 



u 



Figure 15. Z = $ 



u 



tnk ) 



As usual, we are interested in the process given by the eigenvalues of Z as n — > cxd (in our 
setup, k is fixed). Here, with almost the same techniques as in Theorem 16.31 we obtain the 
following result. 
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Theorem 6.5. Almost surely, as n ^ oo, the random matrix ^ <Ei ^{Emk) G ■M.n'^{C) has 
non-zero eigenvalues converging towards 

( \ 



7 



\-t 1-t 



A:2 ' /fc2 



1-t 



V 



fc2_l times 



I 



Proof. Theorem 16.31 is a particular case of this theorem with t = (1 — l/fc) and it has been 
worked out in great detail. Therefore we leave the reader to work out the appropriate modi- 
fications to this case. □ 

The striking fact here is that the largest eigenvalue behaves almost surely like t + {l — t)/k'^. 
The existence of a large eigenvalue was already anticipated by Hayden in [7], Lemma II. 2. 
The lemma below is a slight generalization of Hayden's lemma, following his idea. 

Lemma 6.6. In the above model, the largest eigenvalue is at least t. 

For the sake of being self contained, we give a proof of this fact. Moreover, the proof is 
graphical, using our diagrammatic calculus developed in Section [3l 



Proof. Following [7j, it is enough to prove that 
T^{ZnEn) = (^Bell„, ® 



-'tkn I 



BelL) > t. 



n 



In order to accomplish this, we use the diagram invariance to stretch outside the inner 
parts of the diagram and then we notice that the resulting diagram is of the form {X, X) for 
some X G , see Figure [T6l Introducing the orthogonal projector A = 1^2 —E^ G A^fc2(C), 



Tr(Z„£;„) 



t7ik n 



u 



u 



u* 



tn^k 



X 



u* 



X* 



Figure 16. Re-writing Ti{ZnEn) 



it is obvious that {X,AX) ^ and thus the inequality in Figure [T71 holds. Note that we have 
replaced the identity operator 1^2 connecting X and X* by the maximally entangled state on 

c'^ (g) C". 




Figure 17. Replacing 1^2 by E^ 
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Figure 18. Application of the unitary axioms 



Now we use the unitary axioms on each of the two connected diagrams above and we obtain 
the result in Figure [T51 

Putting all the factors together, we get 

and then, since En is an orthogonal projector, we conclude that the largest eigenvalue of Z„ 
is at least t. □ 

To conclude, let us compare Theorem 16.51 and Lemma 16.61 Independently on the choice 
of t and A;, the value that we obtain almost surely for the largest eigenvalue in Theorem 
16.51 improves strictly the lower bound for the largest eigenvalue obtained in Lemma 16.61 
t + {l- t)/k'^ > t. 

For fixed t, the relative improvement + (1 — t)/k'^)/t becomes small as k becomes big. 
On the other hand, if t ^ the Lemma 16.61 does not bring new information, as the largest 
eigenvalue is always at least /c~^, whereas Theorem 16.51 brings new information. 

So, the relative improvement is optimal for t ~ k~'^. This heuristic study leads us to think 
that it is possible to improve known counterexamples about the various entropy additivity 
conjectures, and this is the object of study of our second paper [3]. 
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